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ABSTRACT 

We compute the leading a s corrections to the two-point correlator of the Oab=2 
operator which controls the B°B° mixing. Using this result within the QCD spectral 
sum rules approach and some phenomenologically reasonable assumptions in the 
parametrization of the spectral function, we conclude that the vacuum saturation 
values Bb — Bb* — 1 are satisfied within 15%. 



CERN-TH.7140/94 
PM 93/16 

KEK Preprint 93-184 
January 1994 



1 On leave from Institute for Nuclear Research of the Russian Academy of Sciences, Moscow, 
Russia 



The violation of the CP symmetry is still one of the most intriguing phenomena 
in particle physics. The standard model (SM) offers a pattern for an explanation 
of this violation through the complex phase of the unitary Cabibbo-Kobayashi- 
Maskawa quark mixing matrix for three generations. The rigorous experimental 
test of the pattern requires knowing the numerical values of some hadronic matrix 
elements at low energies with some non-perturbative QCD methods. 

In this paper we discuss an estimate of the B^-B^ mixing parameter within the 
QCD spectral sum rules (QSSR) approach [1,2] along the lines of Ref. [3], but by 
including the new a s corrections and by taking care on the detailed contributions 
of different B-like states to the spectral function. In so doing, we estimate the two- 
point correlator of the Oab=2 = (&l7^l) 2 operator, where qi are left-handed quark 
fields 

T(x) = (0\TO AB=2 (x)O AB=2 (0)\0). (1) 

The leading term has a trivial expression in the configuration space. Indeed, to 
leading order in a s , it simply reads 

T (x) = 2N 2 (l + -^J 16S'(x,m)S(-x,0)S'(x,m)S(-x,0) 

= 2 (l + i-) tr[S(x, m)S(-x, 0)]tr[S(x, m)S(-x, 0)] = 2 (l + ^) U 5 (x)U 5 (x) 

(2) 

where S(x,m) is the free fermion propagator and iV c stands for the number of quark 
colours. The prime means taking only the part of the propagator that is proportional 
to a 7 matrix. The function n 5 (x) = (0|T , < 7 5 (x)j 5 (0) |0) is the two-point correlator 
associated to the current j 5 = bi^d. 

One can also rewrite Eq. (2) in the form 

T (x) = 2(l + ^U IIU (x)W^x), (3) 

where W v (x) = (0\Tj£{x)j L (0)\0) and f L = b L ^d L , which has the following Lorentz 
decomposition in x-space 

W u {x) = {-d tx d v + g^d 2 )U T (x 2 ) - d»d u U L (x 2 ). 

Formulae (2) and (3) demonstrate an explicit factorization of the correlator (1) in 
the configuration space to leading order in a s . 

The dispersion representation in x-space for any two-point correlator Hj(x) (j = 
T, L, 5) has the form 

ili j {x 2 )= r j (s)D(x,s)ds, (4) 

J S j 

where D( free boson propagator with the "mass" y/s. The spectral functions 

rj, read to leading order in a s : 

r?\s) = N c ^- 2 z{l - z) 2 , r?)( S ) = A^Ifl - z)\2 + z), 
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-fM=m.'W.^2<L-£>_, (5) 

where z = raffs, nib is the o-quark pole mass. 

One can define the spectral function p(s) of the full correlator T(x) in the same 
way as in Eq. (4) and express it, to first order in a s , in terms of the spectral functions 
rj(s) associated to the two- line correlators. Therefore 

p( s ) = J r 1 (s 1 )r 2 (s 2 )®(s;s 1 ,S2)ds 1 ds2, (6) 

where 



$(s; si, s 2 ) = — y (s - si - s 2 ) 2 - 4 Sl s 2 (7) 

is the two-body phase-space factor, and the concrete form of the spectral functions 
Tj(s) entering Eq. (6) depends on the representation chosen for the correlator T(x) 
(as in Eqs. (2,3) for T (x)). The integration region in Eq. (6) is determined by 
the properties of the phase-space factor (7) for corresponding representation of the 
whole correlator; $(s;si,s 2 ) is supposed to be equal to zero within kinematically 
forbidden regions. 

To leading order in l/N c , one can write the correlator T(x) as a product of two 
two-line correlators given in Eq. (4). It is worthwhile to notice that this decomposi- 
tion is gauge-invariant and finite, i.e. it does not require any renormalization that is 
a reflection of the vanishing of the operator Oab=2 anomalous dimension in leading 
order in 1/N C . 

Including a s corrections, the spectral function in the factorization approximation 
reads: 

P fact(s)= p (s)(l + Ap f (s)), (8) 

where po(s) generates the leading-order term and Apf(s) is a (properly normalized) 
factorized correction in the a s order. The colour structure is the following 

Po(s) = Nl (l + A-) p (s), Ap f (s) = C F ^Ap f (s), 

where C F = (iV? — 1)/2N C for the SU (N c ) colour group and the quantities with the 
tilde contain no explicit A^ c dependence. 

The representation of the spectral density Pf ac t{s) through those of two-line 
correlators for the vector-like decomposition in Eq. (3) is: 

Pfact(s) = J ds l ds 2 $(s; s u s 2 ){ + 2s!S 2 ^ r T (s 1 )r T (s 2 ) 

Is 1 \ s 2 

+ ( - sis 2 J (r T (s 1 )r L (s 2 ) + r L (s 1 )r T (s 2 )) + -^-r L (si)r L (s 2 } 

where s i2 = s — si — s 2 . The non-factorizable corrections are of the 1/N C order and 
the full spectral density can be written in the form 

p(s) = p (s) (1 + Ap f (s) + Ap nf (s)) (9a) 
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and 



Ap n f(s) 



C F a s 
N c +1~ 



A*/M = \ (l - £) >P„,M 



(96) 



One of the diagrams contributing to the non-factorizable part of the whole spectral 
density is given in Fig. 1. 

Several comments are in order here. The decomposition (9) is gauge-invariant. 
The same is true for the anomalous dimension of the operator Oab=2- We use four- 
dimensional algebra of Dirac's 7 matrices throughout the computation. This seems 
quite natural because it allows one to make Fierz rearrangements freely, which is 
crucial for establishing the validity of factorization. At the same time this approach 
simply implies the special choice of the renormalization scheme. So, this scheme is 
not the standard MS one and it is not the real dimensional reduction as used in 
[4] either. To the considered order in a s , however, the scheme dependence reduces 
to a certain choice of the normalization parameter fi. In principle, one can fix the 
scheme by direct comparison of our results with the massless limit for corresponding 
correlator [5] or, to put it another way, by comparing the non-logarithmic parts of 
corrections in the massless limit one can obtain the relation between our parameter 
fi and the corresponding fijjg or ^oimRed- Thus, the answer in any desired scheme 
can be easily recovered by considering the massless limit. In the present paper we 
do not dwell on this point. 

We have computed numerically the spectral density p(s) in the first non-leading 
order in a s . We analyse our results, paying special attention to the presentation of 
the entire spectral density p(s) as a sum of factorizable and non-factorizable pieces, 
as in Eq. (9). This decomposition is useful from the theoretical point of view. It 
is also interesting from the 1/N C analysis, and it is quite convenient technically. As 
for the factorizable part of the spectral density (Eq. 8) the analytical expressions 
for two-line spectral functions r TtL (s) are well known to first order in a s , and we use 
them as they have been given in Refs. [6,7]. For the non-factorizable part of the 
spectral density, we have to compute the gluon-exchange diagram both for unequal 
mass lines and for equal mass lines of a two-line correlator. We have done it using the 
REDUCE system of analytical computation [8] and the table of two-loop integrals 
given in [6,7]. After using a Fierz transformation, these non-factorizable diagrams 
can be represented in form analogous to that of the factorizable diagrams, i.e. the 
corresponding analytical expression is given by a product of two traces in Dirac's 
indices. The transformed diagrams are given in Fig. 2. 

For diagrams with a gluon connecting b or rf-quarks in the loop (Fig. 2(a,b)), 
the resulting representation can be rendered into a product of two scalar spectral 
densities of two-line correlators. The relevant part of the spectral density reads: 



r m (s) 



1 



(2(1 - 2z)v + C F 



-Ml-2,)(i + 21n^) 



16tt 2 



s 



7T 



+4(1 



2z)(l 



v)\nz + (l + v)(Azv + (3 + v)(z - 1)) In 



1 -v 



1 + v 
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where e = 2 — D/2, D is the dimensional regularization parameter, v = \Jl — = 
VI - 4i, = 8(Sp(w 2 ) - 5p(-u)) +4(21n(l - u 2 ) - ln(l +u))lnu, u = (1 - 
«)/(! + «), 

JO t 

For a zero mass correlator (rf-quark), the spectral density has the form 
r{s) = -±- 2 s\2 + C F ^ (-4(1 + 2 In ^) - 21^ } . 
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The corresponding representation for diagrams, with different quark flavours 
tied with a gluon in the loop, is of the vector type. The nonfactorizable first order 
corrections to the spectral densities are: 

r L ,T(s)=r^ T (s)+C F ^ T (s), 

71 

r£\s) = z(l - z) 2 {- + 2 In ^4) - 2z(l - ^) 2 ln(l - z) H~^-) 

6 /77^ 1 Z 

+z(l-z) 2 (4z 2 -3z + l-)\n(-^)+4z(l-z) 2 Sp(—-^) + z(l-z)(2z-3)\n(l-z) 

2 1 — z 1 — z 

/, 2 ^ 3 ^ 65 3 17 2 21 2 
z(4z 2 - 72; + -) In z + —z 3 - — z 2 + —z - - , 

3 4 1} ( S ) = (z + 2)(1 - ^) 2 (- + 2 In ^4) - 2(z + 2)(1 - ln(l - z) ln(-^) 

e m% I — z 

+ (1 - z) 2 (2z 3 - 2z 2 + §z + 2) hi(— ^-) + 4(z + 2)(1 - z) 2 ^--^— ) 

+z(l - z)(2z + 1) ln(l - z) + (3z 3 + 2z 2 - \z + 2) In* + ^-z 3 + \z 2 - ^-z + ^. 

The whole computation can be subjected to a powerful test consisting in the 
cancellation of divergences by the one-loop renormalization constant of the operator 
Oab=2, which is known to be Z 0ab=2 = 1 — 3(1 — 1/JV C )^. We checked this 
cancellation explicitly. 

The QCD results for the spectral density are given in Table 1. 

At larger s (say, s > 5m 2 ), where the perturbative QCD result makes sense, one 
can notice that the non-factorizable correction dies out faster than the factorizable 
one, so that the full correction comes mainly from that of the two-line correlator. 
Therefore, at sufficiently large s, the factorization (in the sense of Eqs. (2) and (3)) 
becomes exact and the 'formal' 1/N C suppression of non-factorizable corrections 
becomes numerically valid. 

Using the QSSR approach for estimating B B , we can work with the moments: 

r s th 

M^stn) = / p(s)s- l ds, (10) 

J Ami 



which can be decomposed as: 

Mi( Sth ) = M?(s th ) (l + AM/ (s th ) + AMf(s th )) , 

according to the corresponding decomposition of the spectral density, Eq. (9a). 

We show in Table 2 the strength of the non-factorizable correction for s > 5m 2 h 
for the following set of input parameters: = 175 MeV, m b = 4.6 GeV, and with 

the one-loop expression for the strong coupling constant a s (fj) = 67r/231n(///A^) 
at the point ji = m\>. For given regions of integration in Eq. (10), the moments of the 
factorizable spectral density are practically independent of the power % of the weight 
function s~ l . They change by less than 5% for i=l- 10. The non-factorizable correc- 
tion changes its sign within the integration region for s > 5.5ml, an d its moments 
are more sensitive to the power i. One can notice that the non-factorizable correc- 
tion does not exceed a 15% level with respect to the full factorized spectral density 
at sufficiently large s. The full order-o; s corrections (factorizable+non-factorizable) 
can be quite large for both the spectral density itself and its moments depending 
on the energy s as they can reach a magnitude of 100% with respect to the leading 
term. Hopefully, the non-factorizable corrections, measured in terms of the fully fac- 
torized (lowest order + a s terms) spectral density are, nevertheless, still moderate. 
However, all the factorizable corrections to the correlator (1), no matter how large 
they are, can be absorbed into the calculation of the decay constant fs from the 
two-point correlator with two quark lines, in such a way that the relevant corrections 
to the Bb parameter are only due to the non-factorizable ones. 

Contrary to the QCD part, the estimate of the phenomenological contributions 
to the correlator is much more involved. For the analysis to be performed one might 
make the following phenomenological assumptions: 

- the minimal choice of operator: Oab=2 = 9AB=2d^Bd^B, can give a good descrip- 
tion of the spectral function. 

- the contributions of the BB and B*B* pairs to the spectral function are equal due 
to their approximate degeneracy and to the equal values of their decay constants [9] 
(this feature is fully satisfied in the large mj limit); 

-we do not have exotic contributions due to a singlet mixture of colored states or due 
to some four-quark-like states which may restore the ^-structure of the correlator 
in order to match the proper (l + l/iV c ) construction entering the standard definition 
of the Bb parameter. 

Within these phenomenological assumptions, we can conclude that the effect of 
inclusion of the perturbative corrections into the QCD part of the correlator on 
the Bb parameter is reasonably small and the change is not larger than 15%. The 
absolute value of Bb however cannot be unambiguously established using the pertur- 
bative part of the correlator only. Assuming that non-perturbative non-factorizable 
corrections are small and the factorizable corrections are properly taken into account 
through J'b we find after correcting for exotic contributions in theoretical part of 
the correlator 

B B ^B b * (1.00 ±0.15) 
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where we have not explicitly written the small effect due to the anomalous dimension 
of the operator Oab=2- 

An improvement of our result needs an explicit quantitative check of the previous 
phenomenological assumptions. Most probably, the QSSR approach based on the 
three-point function [10] may be more appropriate, as it does not need some of the 
assumptions which we have used for the two-point function. Unfortunately, one has 
to face, in this highly involved computation of the non-recursive three-loop 

non-factorizable diagrams in order to obtain the a s corrections. 
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s/ml 


A Pf 




A Pn/ /(l + A P/ ) 


5.5 


1.03 


0.02 


0.01 


6.0 


0.95 


0.21 


0.11 


6.5 


0.89 


0.29 


0.15 


7.0 


0.84 


0.32 


0.17 



Table 1: Normalized spectral densities 
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Sth/m'i 


AM f 


AM nf 


AM nf /(l + AM f ) 





5.5 


1.07 


-0.16 


-0.08 




6.0 


0.99 


0.11 


0.06 




6.5 


0.93 


0.23 


0.12 




7.0 


0.88 


0.29 


0.15 


5 


5.5 


1.08 


-0.21 


-0.10 




6.0 


1.00 


0.08 


0.04 




6.5 


0.94 


0.20 


0.11 




7.0 


0.89 


0.27 


0.14 


10 
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1.09 


-0.27 


-0.13 




6.0 


1.01 


0.02 


0.01 




6.5 


0.96 


0.16 


0.08 




7.0 


0.91 


0.23 


0.12 



Table 2: Normalized moments of spectral densities 
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